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Abstract 

Context. Recent calculations of pulsation modes in rapidly rotating polytropic models and models based on the Self- 
Consistent Field method (MacGregor et al. 2007) have shown that the frequency spectrum of low degree pulsation 
modes can be described by an empirical formula similar to Tassoul's asymptotic formula (Tassoul 1980), provided that 
the underlying rotation profile is not too differential (Lignieres & Georgeot 2008, Reese et al. 2009). 
Aims. Given the simplicity of this asymptotic formula, we investigate whether it can provide a means by which to 
identify pulsation modes in rapidly rotating stars. 

Methods. We develop a new mode identification scheme which consists in scanning a multidimensional parameter 
space for the formula coefficients which yield the best-fitting asymptotic spectra. This mode identification scheme is 
then tested on artificial spectra based on the asymptotic formula, on random frequencies and on spectra based on full 
numerical eigenmode calculations for which the mode identification is known beforehand. We also investigate the effects 
of adding random frequencies to mimic the effects of chaotic modes which are also expected to show up in such stars 
(Lignieres & Georgeot 2008). 

Results. In the absence of chaotic modes, it is possible to accurately find a correct mode identification for most of 
the observed frequencies provided these frequencies are sufficiently close to their asymptotic values. The addition of 
random frequencies can very quickly become problematic and hinder correct mode identification. Modifying the mode 
identification scheme to reject the worst fitting modes can bring some improvement but the results still remain poorer 
than in the case without chaotic modes. 
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1. Introduction 

Many stars w ith intermediate or high masses are rapid 
rotators (e.g. iReese etldl 120081 and references therein). 
Rapid rotation causes a number of additional phys- 
ical phenomena which make it much more difficult 
to model the structure and evolution of these stars. 
These include centrifugal deformation, gravity darken- 
ing, baroclinic flows, vario us forms of tur bulence and 
transport phenomena (e.g. iRieuto rd 2006a). Much the- 
oretic al work has gon e into modelling these stars 



., .Mevnet fe Maederl Il997l. iRoxburghT 12004 f2 006. 
Jackson et all 120051. IMacGregor et all l2007t iRieutordl 
2006U lEspinosa Lara fe Ricutord 2007). Naturally, such 



models are subject to uncertainties and therefore require 
observational constraints. Asteroseismology, the study of 
stellar pulsations, is currently the best way to probe the in- 
ternal structure of stars and therefore to constrain such 
models. A number of recent works have therefore fo- 
cused on the effects of rapid rotation, and in particu- 
lar stellar deformation, on stellar pulsations. For acoustic 
modes, th ese include studies based on full eigenmode cal- 
culations (lEspinosa et al.l 12004 lLovekin fc Deupreel [2 008. 



Lovekin et al. 2009L Lignieres et al. 20061 Reese et all 
2006L 120091) and studies ba sed on ray dynamics 
(|Lignieres fe Georgeotl2008l . l2009D . There are also a number 
of works on other types of pulsation modes. On the obser- 
vational side, the CoRoT mission is providing stellar pulsa- 



tion data with unprecedented accuracy. However, in order 
to exploit such data, it is necessary to correctly match the- 
oretically calculated pulsation modes with observed ones. 
This process is known as mode identification. 

Until now, it has been very difficult to identify pulsatio n 
modes in rapidly rotating stars (e.g. iGoupil et al.l [2005). 
This is because mode identification requires a proper under- 
standing of the frequency spectrum of these stars. Such an 
understanding has only been reached recently for a c oustic 
mo des. Using ray dynamics. ILignieres fe Georgeotl (2008) 
and ILignieres fc Georgeotl (|2009l ) recently showed that the 
acoustic spectrum of a rapidly rotating star is a super- 
position of spectra from different classes of modes, each 
with their own frequency organisation. The main classes 
are island, chaotic, whispering gallery modes and modes 
corresponding to a periodic orbit of period 6. Of particular 
interest are the island modes. These are the most visible 
of the regular modes since they are the rotating counter- 
parts to modes with low Z — |m| values. Their fr e quenc y 
organisation has been studied in ILignieres et alj (12006). 
ILignieres fc Georgeotl (|2008l I2009D and lReese et al.l T2008) 
for polytropic models and in lReese et al. ( 2009) f or models 
based on the Self-Consistent Field method ( Jackson et alj 
120051 IMacGregor et al.ll2007l ). A new asymp totic formula, 
similar to Tassoul's formula (jTassoull Il980l). was derived 
involving a new set of quantum numbers based on the ge- 
ometry of these modes. This naturally raises the question as 
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to whether such a formula can be used to identify pulsation 
modes. 

In order to address this question, we develop a new 
mode identification scheme which is described in section 2. 
In section 3, we run an initial series of tests on the mode 
identification scheme using frequencies based on the asymp- 
totic formula, random frequencies and numerical frequen- 
cies based on full 2D eigenmode calculations. This is then 
followed by tests using composite spectra in which random 
frequencies have been added to a set of numerical frequen- 
cies. The final section concludes by discussing the results. 

2. A new mode identification scheme 

Various methods for identifying pulsation modes or de- 
tecting underlying regularities in frequency spectra have 
been invented ove r th e past few years. For instance, 
iBreeer et all (|l999l) and iBreger et ail (|2009t) have worked 
with histograms of frequency differences in order to in- 
terpret pulsation spectra of S Scuti stars. Other similar 
techniques include calculating the autocorrelation func- 
tion or the Fourier transform of the power spectrum to 
identify the large frequency sep aration in solar- like pul- 
sators (e.g. IChaplin et alj [20081 and references therein). 
These procedures yield information on the structure of 
the frequency spectrum rather than a detailed identifica- 
tion for each pulsation mode. Another type of approach 
has consisted in directly comparing the set of observed 
pulsation modes to numerically calculated frequency spec- 
tra from models in a large parameter space. Given the 
computational cost involved in computing each pulsa- 
tion spectrum, various methods have been created in or- 
der to search thro ugh parameter space in a n intel ligent 
way. For instance. | Metcalfe fc Charbonneaul ([2003) and 
ICharpinet et a l. (2005) have used genetic algorithms to find 
best matching mo dels for white dwarfs and sdB stars, and 
iBazot et~a l. ( 2008) have used a Monte Carlo Markov Chain 
(MCMC) approach when studying a Cen A. As opposed to 
other methods, this type of approach yields potential mode 
identifications for each frequency in the observed spectrum. 

The mode identification scheme described here com- 
bines some of the characteristics of the previous methods. 
On the one hand, it would be nice to get more detailed 
information than what is available from histograms of fre- 
quency differences or autocorrelation functions and Fourier 
transforms of power spectra. Furthermore, such methods 
may fail in rapidly rotating stars because the geometric 
term — m£l due to advection of modes by rotation may 
give rise to frequency spacings which are comparable to 
the separation between island modes with consecutive n 
values (which corresponds to half the large frequency sep- 
aration). On the other hand, it is not currently feasible to 
calculate complete pulsation spectra of rapidly rotating 2D 
models at each point in a large multidimensional parameter 
space, even with schemes such as genetic algorithms or the 
MCMC approach. A useful compromise is therefore to use 
an asymptotic formula to calculate approximate frequency 
spectra and to explore the parameter space based on the 
coefficients of this formula in search of best fitting spectra. 

In what fo llows, we us e d the following asymptotic for- 
mula based on lReese et al.l (|2009j) to construct approximate 
frequency spectra: 



where A^, A ; ~, A„, and a are coefficients which depend on 
the stellar structure, £1 the rotation rate and h, I and m 
quantum numbers. These quantum numbers correspond to 
the number of nodes along and parall el to the un derlying 
ray paths (see, for example, Fig. 3 of Reesel l2008f ) and to 
the usual azimuthal order, respectively. In what follows, we 
will refer to h as a radial order although it roughly corre- 
sponds to twice the usual (spherical) radial order. This for- 
mula is an approximation which is valid at low azimuthal 
orders of a more complete formula. Given the computa- 
tional efficiency of calculating a single spectrum using the 
asymptotic formula, a simple scan of the parameter space, 
in which the coefficients and rotation rate are treated as 
the parameters, is performed rather than applying a more 
sophisticated search algorithm. 

A number of preliminary choices are made before apply- 
ing this method. Suitable ranges of h, I and m values must 
be selected. These will determine which frequencies are cal- 
culated in the asymptotic spectra. Also, bounds must be set 
for the parameter space. Besides these choices, the set of 
observed frequencies is assumed to be sorted in ascending 
order as this is needed when matching these frequencies 
with those from the asymptotic spectra. The following pro- 
cedure is then applied to each point in the parameter space: 



1 . An artificial spectrum based on the asymptotic formula 
is created using the values of the different parameters 
and the ranges chosen for n, I and m. 

2. The artificial spectrum is sorted in ascending order us- 
ing a heapsort method. 

3. The observed frequencies are matched to neighbouring 
asymptotic frequencies using dichotomy. When two or 
more observed frequencies are nearest to the same ar- 
tificial frequency, only the first one is matched to the 
frequency, the others being matched to following fre- 
quencies. In some cases this can produce sub-optimal 
solutions. 

4. This match between the artificial and observed frequen- 
cies yields a mode identification. Based on this mode 
identification, the formula coefficients are then recalcu- 
lated through a least-squares minimisation of the stan- 
dard deviation between the artificial frequencies and the 
observed ones. 



Figure [T] illustrates these different steps. 
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Figure 1. Schematic diagram which shows how one asymp- 
totic spectrum can be compared with a set of observed fre- 
quencies to produce a plausible mode identification. 
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The search domain is in terms of the following param- 
eters rather than the original coefficients from the asymp- 
totic formula: 



_Aj- a_ 
A- ' A- ' A- ' A- ' 



(2) 



The advantage of working with these parameters is that the 
artificial spectrum does not need to be sorted again when 
the value of A n or is modified. 

It is important to note that the synthetic spectra are 
periodic with a period equal to Aft, except for the ends 
where frequencies will be missing due to the cutoff in h 
values. As a result of this, it is only necessary to test the 
parameter ttt within a unit interval. Choosing the wrong 
interval will offset all of the n values by a fixed amount. 
Since this is probably one of the least robust parameters, it 
is likely that the mode identification scheme will only yield 
relative rather than absolute n values for real stars. 



3. Results 

3.1. Initial tests 

A number of artificial spectra in which the identification 
of the frequencies is known beforehand were used to test 
the mode identification scheme. These spectra included 
artificial frequencies which follow exactly the asymptotic 
formula, random frequencies, and frequencies of low and 
high order modes of a 25 Mq ZAMS model rotating uni- 
formly at 60% of the critical r otation rate (for d etails on 
the eigenmode calculations, see I Reese et al.|[2009t ). Results 
from these tests are summarised in Table [T] In cases 1, 2, 4 
and 5, 50 frequencies were randomly selected with I values 
equal to or 1 and m values between —3 and 3. The h 
values are given in the fourth column of Table [1] For cases 
1, 2 and 4, this corresponds to 50 out of 84 possible modes, 
whereas there are 224 possible modes for the last case. Case 
3 corresponds to a set of 50 frequencies with random values 
in the same frequency range as case 2. As such, these fre- 
quencies have no corresponding identifications. The third 
column gives the standard deviation Su> between the fre- 
quencies and their asymptotic approximations, normalised 
by Aft. This deviation is defined as follows: 



Slo 



\ 



1 N 



asymp. \ 2 



(3) 



i=i 



where N is the number of observed modes, the "ob- 
served" frequencies and o^ symp ' their asymptotic approxi- 
mations, as based on Eq. {2]). T he coefficients in Eq. ^ can 
be calculated in several ways. iLignieres fc Ceorgeotl feOOS) 
give theoretical formulas for A n and Ar based on travel- 
time integrals of underlying ray paths. However, similar 
formulas for the remaining coefficients are not currently 
available. A more heuristic approach is to calculate a set 
of numerical frequencies and find the corresponding coeffi- 
cients using a least-squares fit. This then raises the issue as 
to which frequencies are to be included in the set. In the 
current context, the most logical choice is the set of "ob- 
served" frequencies specific to each of the cases. Indeed, the 
mode identification scheme can only find the asymptotic 
coefficients based on the frequencies which are available. 
Furthermore, choosing these frequencies yields the lowest 



value for 5uj. This implies, however, that the asymptotic 
coefficients will be different for cases 2, 4 and 5 (as can be 
seen in Table [2]) in spite of the fact that these correspond 
to the same model. 

The mode identification scheme was then applied using 
the same I and m ranges. The range on n used by the scheme 
was chosen to be larger than the range used to generate the 
observations (except for the first case), as can be seen in the 
fifth column of Table [1] The bounds of the parameter space 
are given in Table [5] along with the values of the parameters 
corresponding to the exact solution. The parameter space 
was discretised using 50 uniformly distributed grid points 
in each direction, thus yielding a total of 3.125 x 10 s combi- 
nations, except for case 5 where 100 uniformly distributed 
points were used in each direction. The average computa- 
tional time was approximately 1 hour on a single PC for 
cases 1-4, and around 20 hours for case 5. 



Table 2. Bounds to parameter space and solutions. 



Case(s) 


(in /iHz) 


A f Am 

Aft Aft 


tl 

Aft 


a 

Al 


Bounds to parameter space 


1 


11.5-19.2 


0.5-0.9 0.01-0.05 


0.6-1.0 


2.5-3.5 


2,3 


9.6-19.2 


0.5-0.9 0.00-0.05 


0.8-1.2 


2.8-3.8 


4,5 


9.6-19.2 


0.5-0.9 0.00-0.05 


0.8-1.2 


1.0-2.0 


Solutions 


1* 


17.26 


0.660 0.0288 


0.827 


2.92 


2 


15.25 


0.767 0.0205 


0.955 


3.35 


4 


16.01 


0.755 0.0076 


0.919 


1.64 


5 


16.01 


0.758 0.0074 


0.921 


1.65 



Bounds of the parameter space used in the mode identi- 
fication scheme and corresponding exact solutions. These 
bounds were chosen so as to include the solutions. 
*These parameters correspond to a polytropic model with 
a polytropic index of 3 and the same mass and equatorial 
radius as the model used in cases 2, 4 and 5 (i.e. M — 25Mq 



and i? eq = 7.467? ). 



Column 6 of Table Q] gives the success rate in correctly 
identifying modes, averaged over the 100 best solutions 
given by the identification scheme. The last column gives 
the lowest standard deviation between the fitted frequen- 
cies and observations, normalised by the A n of the output 
solutions. 

As could be expected, using frequencies which follow 
exactly the asymptotic formula yields the best results. 
Although the exact solution is within the bounds of the 
parameter space tested, it was not actually found - the 
nearest grid points yielded mode identifications which were 
slightly different than the original. 

Using low order numerical frequencies does not yield 
good results, as can be seen from case 2. The reason for 
this failure seems straightforward. At low radial orders, the 
deviations caused by avoided crossings between the frequen- 
cies and their asymptotic values can be substantial. As a re- 
sult, the mode identification scheme found erroneous iden- 
tifications which actually led to closer fits to the numerical 
frequencies than the original identifications. This point is 
further confirmed by case 3, in which the mode identifica- 
tion scheme is able to reproduce a set of random frequencies 
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Table 1. Description of the "observed" frequencies and parameters/results from the mode identification scheme. 



. . . . "Observed" frequencies Mode identification 



Case 










Average success 






Type 


A~ 


n 


n 


A~ 




1 


Asymptotic 





10-25 


10-25 


77.4 % 


3.7 x 10" 


3 


2 


Numerical 


2.4 x 1CT 2 


15-20 


10-25 


17.5 % 


1.2 x 10~ 


-2 


3 


Random 






10-25 




2.1 x 10" 


-2 


4 


Numerical 


2.8 x 10" 3 


45-50 


40-55 


65.9 % 


5.4 x 10~ 


-3 


5 


Numerical 


7.9 x 10" 3 


40-55 


35-60 


73.9 % 


7.3 x 10' 


-3 



Characteristics of the "observed" frequency spectra (columns 2-4), of the n values used in the mode identification scheme 
(column 5) and of the results (columns 6-7). The observed frequencies from cases 2, 4 and 5 come from a 25 Mq model, 
uniformly rotating at 60% of the critical rotation rate. 



with no underlying regularities with a standard deviation 
equal to 0.02A,~ l . This is marginally worse than the stan- 
dard deviation obtained in case 2. 

Going to higher radial orders substantially reduces the 
deviations between the numerical frequencies and their 
asymptotic values. This can be seen by comparing the stan- 
dard deviations of cases 2 and 4 (see 3 rd column of Table [T]) 
which shows a tenfold decrease at higher radial orders. This 
leads to good results when applying the mode identification 
scheme. Figure [2] shows a series often plots arranged in tri- 
angular form which give an idea of the accuracy of the fits in 
parameter space. Each plot corresponds to different pairs of 
parameters, which we shall generically denote as p\ and p2- 
The plots show the standard deviation, Suj, as a function of 
Pi and P2, the remaining variables being optimised over the 
parameter space. The colour bar on the righthand side indi- 
cates the meaning of the different colour levels. Since these 
plots are based on the parameters after the least-squares 
adjustment, the different positions in each plots are actu- 
ally bins in which only the best solution is retained. Using 
the adjusted parameters also means that some regions will 
be avoided, and these are indicated in yellow. 

Superimposed on these plots are the positions of the 
100 best solutions, as represented by the red crosses. As 
can be expected, these crosses are concentrated in the dark 
regions which correspond to the best fits. The white di- 
amonds show the exact solution. The plot in the upper 
right corner, besides the colour bar, indicates the number 
of modes correctly identified. As can be seen from these 
plots, the correct solution has a basin around it which at- 
tracts most of the 100 best-fitting solutions. Some other 
secondary basins around other solutions also appear and 
attract a few of the best-fitting solutions. 

The last case in Table[T]tests the effects of using a larger- 
range of n values and therefore a sparser set of observed fre- 
quencies. In order to obtain good results, it was necessary 
to use a finer grid, hence the reason why 100 rather than 50 
points were used in each direction. This need for a higher 
resolution indicates that the basin around the exact solu- 
tion is smaller, probably as a result of the larger range of n 
values. Once a sufficient grid resolution is chosen, then the 
mode identification scheme can yield very good results in 
spite of the sparseness of the set. 

Another interesting test consists in changing the num- 
ber of observed modes. This has been done for case 4 of 
Table [TJ Two different measurements of the success rate 
are shown in Fig. [3j The solid line shows the average suc- 
cess rate of the 100 best-fitting solutions whereas the dot- 



ted line shows that of the best solution. As the number of 
observed frequencies decreases, many additional local min- 
ima appear in parameter space due to the lower number 
of constraints. These cause the average success rate to de- 
crease, especially when there are fewer than 30 modes, by 
attracting an increasing number of the 100 best-fitting so- 
lutions away from the true solution. The dotted line, on the 
other hand, shows that in most cases, the best-fitting solu- 
tion still remains the true solution. However, in some cases, 
one or several other local minima with completely different 
mode identifications produce solutions which fit the obser- 
vations even better than the true solution, thereby causing 
a dramatic drop in the accuracy of the best solution. This 
is illustrated for 11 and 16 observed modes in Fig. [31 As a 
result, the best-fitting solution must be used with caution. 



ioo ft^: 



^ 80 r ; 1 ; 




10 20 30 40 50 

Number of observed modes 



Figure 3. The average success rate of the 100 best-fitting 
solutions (solid line) and the success rate of the best solu- 
tion (dotted line). 



3.2. Tests with additional random frequencies 

As was shown in lLignieres fc Georgeotj (|2009T ) based on geo- 
metric visibility calculations, chaotic modes are likely to be 
visible in the pulsation spectrum of rapidly rotating stars. 
The frequencies of these modes do not follow an asymp- 
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Figure 2. Plots which show the stan- 
dard deviation, Slu, in parameter space 
for case 4 of Table [TJ The ten plots ar- 
ranged in triangular form correspond to 
different pairs of parameters. For each 
pair, the plot shows 5ui as a function of 
the two parameters, the remaining pa- 
rameters having been optimised over the 
search domain. Superimposed on these 
plots are the 100 best solutions (as shown 
by the red crosses) and the exact solution 
(indicated by the white diamond). The 
plot in the upper right corner shows the 
number of correctly identified pulsation 
modes for the 100 best solutions. 



10 12 14 16 18 0.5 0.6 0.7 0.8 0.00 0.01 0.02 0.03 0.04 0.8 0.9 1.0 1.1 

A„ (in ,tiHz) Aj/A,, A m /A n n/A„ 



totic formula but rather a statistical distribution. These 
will naturally make mode identification more difficult as it 
is not possible to know a priori which modes are regular 
and which ones are chaotic in an observed frequency spec- 
trum. In order to mimic the presence of chaotic modes, we 
have done a number of tests in which random frequencies 
were added to the frequencies from cases 2 and 4 of Table Q] 
Figure [4] shows two plots with the average success rate 
as a function of the number of additional random frequen- 
cies, one for case 2 (left panel) and the other for case 4 
(right panel). Although low order modes yielded worse re- 
sults without random frequencies, they seem to be less af- 
fected by the presence of random frequencies than high or- 
der modes. The reason why the results for high order modes 
are so poor seems to be that random frequencies shift the 
basin of best-fitting solutions away from the true solution. 
Even a small deviation between the two can be sufficient to 
throw off the identification of modes. 

One way of trying to deal with random frequencies is to 
add another step in the mode identification procedure. At 
each point in parameter space, after matching the artificial 
frequency spectrum with the observed frequencies (step 3 of 
the procedure described in Section 2), one can remove the 
N cut worst modes, where A cut is a fixed number. Figure 
illustrates the mode identification procedure with this extra 
step. 

Figure [6] shows the effects of removing A cut = 20 fre- 
quencies on the average success rate as a function of the 
number of additional random frequencies. The frequency 
sets are the same as those used in the right panel of Fig. [4] 
The success rate has been calculated only using the modes 
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Figure 5. Same diagram as Fig. [T] except for an additional 
step in which the A cut worse frequencies are removed. 



which have been retained: rejecting non-random frequencies 
has not been penalised. As can be seen, having a compara- 
ble number of random frequencies to N cut yields the best 
results. The reason why the results are poor when there are 
no random frequencies is because additional solutions with 
lower standard deviations than the original solution are ap- 
pearing. This is similar to what happened when the number 
of observed modes was reduced (see Fig. [3]), although this 
time the problem is worse because the mode identification 
is carefully choosing the subset of frequencies so as to re- 
duce the standard deviation. When random frequencies are 
added, these tend to be rejected rather than the true fre- 
quencies, thereby leading to more favourable results. Even 
in the best of situations, identifying modes remains difficult 
when random frequencies are present. 
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Figure 4. Average success rate as a function of the number of additional random frequencies. The left plot corresponds 
to low order modes (case 2 of Table [J) and the right plot to high orders (case 4) . 




10 20 30 40 

Number of random frequencies 



Figure 6. Average success rate as a function of the number 
of additional random frequencies. The frequency sets are 
the same as those used in the right panel of Fig. [4] This 
time, the A cut = 20 worst frequencies have been rejected 
at each point in parameter space. The success rate has been 
calculated based on the modes which have been retained: 
rejecting non-random frequencies has not been penalised. 



4. Conclusion 

In this paper we have created and tested a new method 
for identifying pulsation modes in rapidly rotating stars, 
based on an asymptotic formula for island modes, the ro- 
tatin g counterparts to acoustic modes with low I — \ m\ val- 
ues (jLignieres fc Georgeotl [20081 iReese et all 120091) . This 
method consists in scanning a parameter space in search 
of the formula coefficients which lead to the best agree- 
ment with observations. Results show that it is possible to 
correctly identify pulsation modes provided they are suf- 
ficiently numerous and their frequencies close enough to 
their asymptotic values. Such a situation occurs for 30 or 



more high order modes, typically in the h = 40 — 50 range. 
Bearing in mind that n is roughly twice the spherical ra- 
dial order, this range is similar to that of the solar p-modes. 
This method may apply to S Scuti stars thanks to their nu- 
merous pulsation frequencies, but these will have to be of 
sufficiently high radial order. 

However, when even a few random frequencies are added 
in order to mimic the presence of chaotic modes, results 
can become quite poor. Adding an extra step, in which the 
worst modes are removed, can bring a small improvement 
but this still remains insufficient to enable a reliable mode 
i denti fication. Recent calculations bv lLignieres &: Georgeotl 
(|2009l ) suggest that chaotic modes may be as much as 5 
times as numerous as island modes for a star at 60% of 
the critical rotation rate, which is far more than what is 
considered here. One way to try to deal with this problem 
is to start with stars at lower rotation rates, where a smaller 
fraction of the modes are chaotic. 

One of the weaknesses of this method is the large num- 
ber of input parameters. These include the ranges on the 
different quantum numbers and the bounds on parame- 
ter space. In order to set these parameters correctly, prior 
knowledge of the star through other observations will be 
needed as well as an understanding of the relationships be- 
tween the coefficients in the asymptotic formula and fun- 
damental stellar parameters such as the mass, age and ro- 
tation rate. Such an understanding can only be reached 
through a systematic study of the frequency spectra of a 
grid of stellar models. 

Potential improvements in this method are as follows. 
Rather than doing a simple scan of the parameter space, us- 
ing a more sophisticated search method like a genetic algo- 
rithm (e.o. lMetcalfe fc Charbonneaull2003llCharpinet et all 
120051) or an MCMC method (e.g. iBazot et al.l l2008D would 
enable a more detailed investigation of regions of interest 
while spending less time in other less important parts of 
the parameter space. This would be beneficial to case 4 of 
Table [1] as the solutions that were found were not as op- 
timal as the true solution, and case 5 because the basin 
around the true solution is smaller. Also, correlating mode 
visibilities with observed mode amplitudes could provide 
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supplementary constraints. This, of course, is only viable if 
a simple way to estimate mode visibilities exists. It is then 
hoped that the results from this method would provide a 
useful starting point for comparing full 2D calculations to 
observations. 
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